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ABSTRACT 

We examine the local stability of galactic discs against axisymmetric density perturba- 
tions with special attention to the different dynamics of the stellar and gaseous components. In 
particular the discs of the MiUcy Way and of NGC 6946 are studied. The Milky Way is shown 
to be stable, whereas the inner parts of NGC 6946, a typical Sc galaxy from the Kennicutt 
(1989) sample, are dynamically unstable. The ensuing dynamical evolution of the composite 
disc is studied by numerical simulations. The evolution is so fierce that the stellar disc heats up 
dynamically on a short time scale to such a degree, which seems to contradict the morpholog- 
ical appearance of the galaxy. The star formation rate required to cool the disc dynamically is 
estimated. Even if the star formation rate in NGC 6946 is at present high enough to meet this 
requirement, it is argued that the discs of Sc galaxies cannot sustain such a high star formation 
rate for longer periods. 
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Q 1 INTRODUCTION 

H 

It has been long suspected that some Sc galaxies are so gas rich 
, that their gaseous discs reach the threshold of dynamical instability 

IL" • (Quirk 1972). In an influential study Kennicutt (1989) has shown 
. 5^ for a considerable sample of Sc galaxies by careful analysis of the 
distribution of atomic and molecular hydrogen in each galaxy that 

5_j in the inner parts of the gaseous discs of the galaxies the Toomre 
(1964) stability condition is indeed violated. Kennicutt argues fur- 
ther that, since the galactocentric distances of threshold of dynami- 
cal instability coincide closely with the outer boundaries of the HII 
region discs of the galaxies, dynamical instabilities have led to the 
enhanced massive star formation rate in the inner parts of Sc galax- 
ies. There are, however, counter examples in his sample, notably 
M33 and NGC 2403, which do not reach the threshold level. In ad- 
dition Ferguson et al. (1994) have shown that HII regions can be 
also found in the outer - stable - parts of the discs of galaxies. 

Even though it is intuitively plausible that dynamical instabil- 
ity leads eventually to enhanced star formation as observed in the 
Sc galaxies, we wish to point out that the existence of an unstable 
gas disc has grave consequences for the dynamics not only of the 
gaseous disc but for the stellar disc as well. For this purpose we 
examine in section 2 stability criteria for composite stellar and gas 
discs and derive a dispersion relation in order to estimate the time 
scale on which the instabilities develop. In section 3 we apply the 
stability criterion to the discs of the Milky Way and NGC 6946, a 
typical representative from the Kennicutt (1989) sample. The disc 
of NGC 6946 is shown to be dynamically unstable and we study 
its fierce dynamical evolution by numerical simulations in section 



4. The implications of these simulations are discussed in the final 
section. 



2 STABILITY CRITERION 

The dynamical stability of the composite gaseous and stellar disc 
is examined, which gives a local stability criterion. Spiral instabil- 
ities and the bar instability, which are less violent, are not consid- 
ered in this section. In previous stability studies of multicomponent 
galactic disc models the stellar disc has usually been described by 
leans equations (Biermann 1975, log & Solomon 1984a,b, Bertin 
& Romeo 1988, Elmegreen 1995, Jog 1996). Since the shape of the 
velocity distribution of the stars governs to some degree the stabil- 
ity of the stellar disc we prefer a full stellar-dynamical treatment. 

It has become customary to describe the stellar velocity distri- 
bution by a Schwarzschild distribution. Wielen and Fuchs (1983), 
however, have pointed out that a velocity distribution with an ex- 
ponential shape appears to be more realistic, at least for stars in the 
Milky Way, because the velocity distribution is made up of many 
generations of stars with velocity dispersions varying according to 
their age. First attempts to model such a distribution in the context 
of stability studies have been made by Morozov (1981). 

In Fig. |l| a one-dimensional projection of a distribution func- 
tion of the form 



(1) 



2-Kauav \ \(Tu J \<Jv J 

is shown in comparison with the observed velocity distribution of 
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the McCormick stars - a kinematically unbiased sample - in the 
Third Catalogue of Nearby Stars (Gliese and JahreiB 1991). In 
equation (1) U and V are the radial and circumferential veloc- 
ity components of the stars, respectively, and au and av are the 
corresponding second moments of the velocity distribution which 
we call loosely in the following velocity dispersions. The ratio of 
the velocity dispersions is given by the epicyclic ratio av / au = 
—2B/Hi, with hi the epicyclic frequency and B Oort's constant. As 
can be seen in Fig. |l] there are a lot of low velocity stars, which is 
well modelled by the distribution function (1). 

In his original study of the stability of stellar discs against ax- 
isymmetric perturbations Toomre (1964) assumed a Schwarzschild 
distribution for the stellar velocities. Graham (1965) and indepen- 
dently Toomre (cited in Graham (1965)) generalised the results to 
arbitrary distribution functions depending on the radial action inte- 
gral U"^ + -^g^V"^ (see also the appendix). The distribution func- 
tion (1) is of this type. Graham (1965) did not consider explicitly 
the case of an exponential distribution. But Toomre has treated this 
case in his unpublished material, which he kindly made available 
to us. The main result is that in case of an exponential distribution 
function the standard Toomre stability parameter is slightly modi- 
fied as 



15 - 



(2) 



3.944GEo 

where Eo denotes the surface density of the disc and G is the con- 
stant of gravitation. The numerical factor 3.944 replaces a factor 
3.36 in the case of a Schwarzschild distribution or a factor n in the 
case of an ideal gas. This indicates that a stellar disc with an ex- 
ponential distribution is slightly more unstable than a disc with a 
Schwarzschild distribution, which has fewer low velocity stars. 

We treat the interstellar gas as an isothermal gas, which will 
be sufficient to find the transition of the composite disc to gravita- 
tional instability (Cowie 1981). Graham (1967) has indicated how 
to extend his or Toomre's work to multicomponent disc models. 
Unfortunately, neither he nor Toomre have treated the case of an 
isothermal gas disc inbedded in a stellar disc with an exponen- 
tial velocity distribution, which we consider here. In the appendix 
we discuss the stability of such a disc against axisymmetric per- 
turbations adapting the analyses of Toomre and Graham to the 
present problem. The domain of neutrally stable perturbations in 
parameter space, i.e. the stability parameters of the gas and stel- 
lar discs, respectively, and the gas-to-stellar surface density ratio, 
indicates then which discs are dynamically stable. For the cases, 
where the stability condition is violated, we derive a dispersion re- 
lation, which allows the determination of the wave length and ex- 
ponential rise time of the fastest growing perturbations. 



3 STABILITY OF GALACTIC DISCS 

The stability criterion developed in section 2 is applied to the discs 
of the Milky Way and NGC 6946. 



3.1 Milky Way 

Table 1 summarises the parameters of the Galactic disc which we 
adopt. 

The resulting stability parameters are Q, = 2.4 and Qg = 2.2, 
respectively. On the other hand, using again the parameters given 
in Table 1, the curve of neutrally stable perturbations derived from 
equation (A 19) may be calculated as shown in Fig. 2. The minimal 
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Figure 1. Distribution of- radial - U velocities of 317 McCormick stars in 
the solar neighbourhood. The observed distribution is approximated by the 
projection of an exponential distribution with a radial velocity dispersion of 
48 km/s. 



Table 1. Local Pai ameters of the Galactic Disc 



(TU- 



36.6 

48 
50 
7.5 
6 



km/s/kpc° 
km/s*" 

km/s 



lAU 1984 galactic constants and assuming a 
flat rotation curve 

from solar neighbourhood data (cf. Fig. |l]) 
Kuijken and Gilmore 1989a,b,c 
Dame 1993, multiplied by a factor of 1.4 to 
account for heavy elements 



value of the stability parameter Qg required to stabilise the disc is 
Qg.min = 1.1 implying that the Galactic disc is stable in the solar 
neighbourhood. 

The stability of the inner and outer parts of the Galaxy may be 
analysed in the same way. The surface density and - squared - ve- 
locity dispersion of the stars can be extrapolated by an exponential 
law with a scale length of 4.4 kpc (Lewis & Freeman 1989). The 
surface density of the interstellar gas is taken from Dame (1993) 
and the turbulent velocity dispersion is assumed to be constant. The 
resulting stability parameters are listed in Table 2, implying stabil- 
ity of the Galactic disc against axisymmetric perturbations in the 
inner and outer parts as well. Similar conclusions were drawn by 
Elmegreen (1995) and Jog (1996) on the basis of their models. 

Dynamical stability of the disc is consistent with the observa- 
tion that there is no enhanced massive star formation such as in the 
Sc galaxies of the Kennicutt (1989) sample in the Milky Way. As 
was stressed previously by Ferguson et al. (1994) one finds, on the 
other hand, still star formation at moderate rates in galactic discs, 
which are dynamically stable in the sense discussed here. Appar- 
ently the complex structure of the interstellar medium, which is not 
modelled here, still allows localised Jeans collapse. It appears that 
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Table 2. Stability Parameters of the Galactic Disc 



R = 4.5 kpc R© = 8.5 kpc R = 12 kpc 





2.9 


2.4 


2.6 


Qg 


2.1 


2.2 


1.7 


Qg,min 


1.1 


1.1 


1.3 



Table 3. Stability Parameters during the Evolution of the Galactic Disc 
(R = R0) 





Q* 


Qg 


Qg,min. 


10^° yrs 


2.4 


2.2 


1.1 


5 109yrs 


1.7 


2.2 


1.4 


2.5- 10^ yrs 


1.2 


2.2 


5.5 


1.25- lO'' yrs 


0.9 


2.2 


> 20 



dynamical instability is not a necessary condition for star formation 
but greatly enhances the star formation rate. 

The radial scale length of the Milky Way is not well con- 
strained and, as can be seen in the synopsis by Sackett (1997), may 
be for the main populations of the stars in the Galactic disc as small 
as 2.5 kpc. This would make the inner parts of the Galactic stellar 
disc more unstable. However, the effect is not significant. Assum- 
ing that the Galactic disc has a constant vertical scale height as 
function of galactocentric radius one would expect the - squared - 
velocity dispersion of the stars to have the same radial gradient as 
the surface density of the stars, which implies Q* = 2.1 at i? = 
4.5 kpc. 

The finite thickness of the disc, wliich has not been considered 
here explicitly, has a further stabilising effect on the disc. 

It is interesting to see whether due to the age dependence of the 
velocity dispersion of the stars the Galactic disc was dynamically 
unstable at earlier stages of its evolution. The velocity dispersion 
varies today as crfj oc r, with t denoting the ages of the stars (Wie- 
len 1977). So, assuming a constant star formation rate, it is to be 
expected that at an earlier epoch tep the velocity dispersion of all 
stars which were born until then was 

tep 

"■y = 7— / cru{r) dT=^ af, {t = tep). (3) 



Table 3 shows the evolution of the stability parameters of the Galac- 
tic disc due to this effect assuming that the other parameters of the 
disc were the same as today. 

Apparently, this indicates instability of the disc during early 
stages of galactic evolution. This would be in line with Kennicutt 
et al.'s (1994) argument that early type spirals had apparently an 
initial burst of star formation. However, we believe, for reasons ex- 
plained in the next sections, phases of dynamical instability extend- 
ing over more than 10^ years unlikely. Probably, the surface density 
of the disc was then lower than today and a lot of material might 
have been accreted later onto the disc (Gunn 1982). If, for exam- 
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Figure 2. Curve of neutrally stable axisymmetric perturbations of a gas 
disc embedded in a stellar disc. Local parameters of the Galactic disc have 
been used. 

pie, the surface density was only one half of its present value and 
the gas density equal to the stellar density, the stability parameters 
were Q, = 4, Qg = 1.3, whereas the minimal value of the parameter 
required to stabiUse the disc was Qg^min = 1.2. 

3.2 NGC6946 

NGC 6946 is a typical representative of the Sc galaxies studied by 
Kennicutt (1989). The absolute mass as well as the radial distribu- 
tion of atomic and molecular hydrogen in NGC 6946 have been de- 
termined by Tacconi & Young (1986). The rotation curve has been 
observed by Carignan et al. (1990), who have also constructed a 
mass model comprising a disc and dark halo component. Unfor- 
tunately, the stellar velocity dispersion is not known, but one may 
derive an estimate from the vertical scale height of the disc. Friese 
et al. (1995) have shown in a statistical flattening analysis of faint 
spiral galaxies in the ESO-Uppsala catalogue that the typical in- 
trinsic ratio of vertical to radial scale lengths of Scd galaxies is of 
the order of zo/h = 0.25 implying zo = lA kpc. This may be used 
with the vertical hydrostatic equilibriimi condition to derive an es- 
timate of the vertical velocity dispersion of the stars, 

aw = T-kGTjozo- (4) 

,7^ is a form factor which is equal to 1 for an isothermal disc 
with a vertical sech? density profile and about 2 for the Milky Way 
disc (Friese et al. 1995). The same value is also adopted for the 
disc of NGC 6946. The radial velocity dispersion is estimated by 
using the same axial ratio of velocity ellipsoid as in the Milky Way, 
cr^/a^ = 1 + K^/4B^. In columns 5, 8, and 9 of Tab. 4 the radial 
variations of the stability parameters, which have been calculated 
including finite thickness corrections, are shown. We note that the 
stability parameter of the stellar disc determined in this way, Q* « 
2, is in good agreement with values determined by Bottema (1993) 
for galaxies for which kinematic data are available. 

There is a distinct drop of the stability parameter of the gas 
disc Qg inside R = 16 kpc, which, as was pointed out by Kennicutt 
(1989), is correlated with the outer boundary of the HII region disc 
of NGC 6946. As can be seen from colimm 9 of Table 4 the sta- 
bility criterion as presented here indicates in accordance with Ken- 
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Table 4. Radial Variation of Stability Parameters of the Disc of NGC 6946 (d = 10 Mpc) 



Qg 



^ g,min 



5 

10 
12 
14 


133 
55 
38 

27 


50 
23 
20 
17 


114 

89 
76 
66 


2.5 
2.2 
2.3 
2.4 


38 
21 
17 
14 


6 
6 
6 

6 


0.6 
0.5 
0.5 
0.5 


0.9 
1.1 
1.1 
1.1 


0.9 
1.4 
1.4 
1.4 


0.7 
1.0 
1.2 
1.5 


16 


19 


15 


52 


2.3 


6 


6 


1.0 


1.0 






18 
20 

kpc 


13 
9 

Mq 


13 
12 

km/ s 
kpc 


44 

39 

km/s 


2.6 
3.0 


5 
5 

pc'-^ 


6 
6 

km/s 


1.1 
1.0 


1.0 
1.0 




kpc 



Molecular hydrogen surface densities have been rescaled to the CO density - to - H2 column density conversion factor given by 
Dame (1993). Atomic and molecular hydrogen surface densities have been multiplied by a factor 1.4 to account for heavy elements. 
Corrected for the finite thickness of the discs, assuming 20,* = 1.4 kpc and 20,9 = 50 pc. 



1^ 




crit,g 



Figure 3. Dispersion relation of exponentially growing axisymmetric per- 
turbations of a composite gas and stellar disc. Parameters of the disc of 
NGC 6946 at galactocentric distance R = 10 kpc have been used. 



nicutt's conclusions that the inner parts of the disc of NGC 6946 
are actually dynamically unstable. This is mainly due to the large 
gas content of NGC 6949, which is quite typical for late type spi- 
ral galaxies (Kennicutt et al. 1994). The exponential growth rates 
as well as the most unstable wavelengths, which may be calculated 
from the combined dispersion relation (A19) are illustrated in Fig. ^ 
and are given in columns 10 and 11 of Table 4. The e-folding rise 
time of the instability is rather short, only 3-10^ yrs at R = 10 kpc 
corresponding to 0.1 epicyclic periods. The unstable wavelengths 
are short compared to the radial extent of the galactic disc. 

Using equations (A12) and (A17) one can show that, even if 



the instability grows at the same rate in the gaseous and the stellar 
discs, the amplitude of the disturbance of the stellar disc is much 
smaller than in the gaseous disc. For the parameters adopted in Ta- 
ble 4 the relative star-to-gas density contrast is typically of the or- 
der of 10~^ so that, even when the perturbations of the gas disc 
become highly non-linear, the stellar disc is only mildly perturbed. 
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Figure 4. Surface density distribution of the stellar (dotted line) and the 
gas (solid line) at time step t=1.2 ■10'^ yrs of the simulation. 
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Figure 5. Rotation curve of the stars (dotted line) and the gas (solid line) 
at time step t=1.2 ■ 10^ yrs of the simulation. 



4 NUMERICAL SIMULATIONS 

The onset of the instability into the non-linear regime can be fol- 
lowed by numerical simulations. Simulations of a dynamically un- 
stable gas disc inbedded in a stellar disc were first carried out by 
Carlberg & Freedman (1985). We have rerun similar simulations 
in order to analyse them in the present context. For this purpose 
we have used a code developed by F. Combes and collaborators 
(cf. Casoli and Combes 1982, Combes and Gerin 1985 for the 
full details). The code implements a two-dimensional stellar disc 
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Table 5. Parameter adopted for the simulation 



Figure 6. Stability parameters of the stellar disc (dotted Hnes) and the gas 
disc (solid lines) at time steps t=6-10^ yrs (upper panel) and t=18 ■ 10^ yrs 
(lower panel) of the simulation. 



(Nstar ~ 38 000) into which interstellar gas clouds are embedded 
(^clouds < 38 000). This composite disc is surrounded by a rigid 
bulge and dark halo potential. The gravitational potential of the disc 
is calculated by a standard particle-mesh scheme and the orbits of 
the stars and the clouds are integrated numerically. The stars inter- 
act only by - softened - gravitational forces, while the gas clouds 
may interact inelastically. This is simulated in an elaborate cloud- 
in-cell scheme, which describes the coalescence and fragmentation 
of the gas clouds. In this way a mass spectrum of the gas clouds is 
established. The various parameters of the collision scheme have 
been adjusted by Combes & Gerin (1985) in such a way that the 
mass spectrum of the gas clouds resembles the mass spectrum of 
molecular cloud complexes in the Milky Way. A finite lifetime of 
4-10^ years due to massive star formation is assumed for the gi- 
ant molecular clouds (GMC) at the high end of the mass spectrum, 
Mgmc > 2 • 10^ M0. After that the clouds are disrupted into small 
fragments expanding initially isotropically at relative speeds of 10 
km/s. 

The softening length of the gravitational forces is 500 pc. Fol- 
lowing Romeo (1994) we have adjusted it in this way, because 
on one hand it is still considerably smaller than the critical wave- 
lengths of the interstellar gas and stellar discs, respectively (see be- 
low). On the other hand, it is of the same order of magnitude of the 
expected vertical scale height of the Galactic disc, so that the sta- 
bilising effect of a finite thickness of the disc is taken into account 
to some extent in our two-dimensional simulations. 

The stars are initially distributed according to a Toomre disc 
(Toomre 1963), 



2tv 



-3/2 



(5) 



where Md denotes the mass of the disc and d is the radial scale 
length. The gas clouds are distributed in an exponential disc. The 
halo and the bulge are modelled by Plummer spheres. The param- 
eters which we have adopted for the various components of the 
model are summarised in Table ^. The begin of the simulation is 
illustrated in Figs. ^ to |^ where the surface densities, the rotation 
curve, and the stability parameters are shown. Fig. ^ shows that 
the the combined disc of stars and interstellar gas clouds is dy- 
namically unstable and resembles in that the inner parts of the Sc 
galaxies in the sample of Kennicutt (1989). The next steps of the 
evolution of the disc are shown in Fig. ^ at multiple intervals of At 



stellar disc mass 


1.1 


IOIOM 


scale length d 


4 


kpc 


truncation radius 


32 


kpc 


gas disc mass Mg 


0.6 


10i°M 


exponential scale length 


16 


kpc 


dark-halo mass 


10 


io^om 


core radius 


24 


kpc 


bulge mass Mi, 


2.4 


10i°M 


scale length 


1.5 


kpc 


mass ratio {M^ + Mg)/Mtot 


0.40 





= 1.2- 10^ years which correspond to about 0.1 epicyclic periods. 
As expected ring-like density perturbations appear immediately in 
the gas disc. The wave lengths of these perturbations are of the or- 
der of the critical wave length of the gas disc, Xc-u.g about 2 kpc 
at R= 10 kpc. The rings fragment into lumps with masses in the 
range of 10* to 10^ Mq. These agglomerates are so heavy that the 
stellar disc responds to them by induced, 'swing-amplified' spiral 
structures. According to the 'X=2' criterion (Toomre 1981) swing- 
amplification is most effective for spirals with a number of 



27ri? 

^crit 



(6) 



spiral arms. Since Xcrit = 10 kpc, the expected number of spi- 
ral arms is about m = 3, which is typically seen in our simulations 
(see for instance time step 14.4 • 10 ). The potential troughs of the 
stellar disc, on the other hand, begin to trap much of the interstellar 
gas and during their further evolution both the stellar disc as well as 
the gas disc undergo rather synchronous, repetitive cycles of swing 
amplified spiral perturbations. 

After 5-10** years the stellar disc gets heated up dynamically 
by the spiral activity so much that hardly any non-axisymmetric 
structure is longer possible in the disc. 

In the gaseous disc, however, there is still a lot of spiral ac- 
tivity. Since the stellar disc has become inactive the critical wave- 
lengths of the spiral structures are much smaller and the number of 
expected spiral has risen to about 20 in accordance with flocculent 
appearance of the disc. The same phenomenon was described by 
Carlberg & Freedman (1985) 



5 DISCUSSION 

The numerical simulations show that the discs of Sc galaxies like 
NGC 6946 are in a highly peculiar dynamical state. The reaction 
of the stellar discs to dynamically unstable gas discs is so fierce 
that they become dynamically hot within less than 10^ years. On 
the other hand, Toomre (1990) has argued emphatically that Sc 
galaxies must have stellar discs, which are dynamically active, be- 
cause otherwise the morphological appearance would be quite dif- 
ferent from what is actually observed. This can be clearly seen 
in Fig. 7, when one compares the frames corresponding of say 
1.44 ■ 10* years and 4.8-10** years with an optical image of the 
galaxy (Sandage & Bedke, 1988). Thus the stellar discs in Sc galax- 
ies must be effectively cooled dynamically by newly formed stars 
on low velocity dispersion orbits. The star formation rate required 
to keep the stellar discs in a steady state can be roughly estimated 
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Figure 7. Dynamical evolution of the stellar and gaseous discs. On the left of each panel 19000 out of 38 000 stars and on the right hand side of each 
panell9 000 of about 38 000 interstellar gas clouds are plotted at consecutive time intervals. The time is indicated on the left in units of lO'^ years. The spatial 
size is indicated at the bottom by a bar of 10 kpc length. 
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Table 6. Radial variation of the stellar disc mass increase per 10^ yrs in 

NGC 6946 



R SFR Sgas S* 






285 


325 


306 


0.9 


1.1 


180 


267 


143 


0.7 


1.2 


71 


219 


95 


0.3 


4.4 


57 


148 


53 


0.4 


6.6 


28 


100 


36 


0.3 


8.8 


11 


67 


27 


0.2 


11 


9 


46 


15 


0.2 


13.2 




31 


11 


0.1 


kpc 




Mq 


Mq 





as follows. After 5 ■ 10 years the stellar disc has heated up so much 
that it tunes out of spiral activity and the - squared - stability pa- 
rameter has risen to 

2 2 

^* " (3.36GE*)2 

Newly formed stars will lower this to 

,2 _ «:^(E*crJ + 5*S*crJ,i) 

^* " (3.36G)2(l + 5.)'Sr 

where we assume that the surface density of the stellar disc has 
risen due to star formation to = Et(l + 5*) and that the - 
squared - velocity dispersion of the mixture of older and younger 
stars can be approximately estimated by an mass weighted average 
of the velocity dispersions of the components. Combining equation 
and (|) gives 

(l + ..f^f(l + ^.^)-f, (9) 

if the velocity dispersion of the newly bom stars, cr^^^i, is much 
smaller than the average velocity dispersion of the stars. In our sim- 
ulations the stability parameter Q* rose from about initially 2 to 
2.5 at t = 4.8 ■ 10^ yrs. According to equation (9) about 40 % of 
the mass of the stellar disc is required per lO'"* years in the form of 
newly bom stars to keep the stellar disc in a steady dynamical state. 
Obviously this material must be provided by the gaseous disc. It is 
interesting to compare this estimate of the gas consumption rate 
with the actual gas consumption rate in NGC 6946. This can be 
deduced from the star formation rate which in tum can be quanti- 
tatively estimated from the Ha flux using the relation of Kennicutt 
(1983). In Table 6 we give the radial distribution of the star for- 
mation rate determined from the extinction corrected Hq, surface 
emissivity of NGC 6946 (Devereux & Young 1993). The surface 
densities of the gaseous and stellar discs shown in Table 6 are taken 
form Table 4. In the last column of Table 6 we give the relative in- 
crease of the stellar surface density per 10^ years, 5*, estimated 
from the previous parameters. As can be seen from Table 6 the ac- 
tually observed star formation rate and thus the disc mass increase 
is at present, if compared with the theoretical estimate found above, 
high enough to keep the stellar disc dynamically cool. In this aspect 
the disc seems to be self regulating. However, the gas consumption 
rate is very large. As is shown in Table 6 star formation seems at 
present to consume nearly the entire gas disc within 10^ yrs. This 
implies either a very high gas accretion rate, which at this magni- 
tude seems unlikely to us, or that the discs of the Sc galaxies of the 



NGC 6946 type will soon switch over to more quiescent dynamical 
states like in M33 or NGC 2403. In these galaxies the gaseous discs 
do not reach the threshold of dynamical instability. But even if the 
stability parameter Qg is slightly larger than 1, this would mean 
still quite a lot of spiral activity in the disc. This might account for 
the still comparatively high star formation rates observed in these 
galaxies (Kennicutt et al. 1994), but leads also to considerable dy- 
namical heating (Carlberg & Sellwood 1984). The same comment 
applies to the low surface brightness galaxies discussed by Mihos 
etal. (1997). 

A further interesting example is the Sb galaxy NGC 7331. Its 
radial distribution of atomic and molecular hydrogen is given by 
Young & Scoville (1982). After rescaling the molecular hydrogen 
densities to a modern CO intensity-to-H2 column density ratio 
(Dame 1993) it can be shown that the gaseous disc of NGC 7331 
does not reach the threshold of dynamical instability, Qg>i. In 
accordance with this result the star formation rate seems to be 
in NGC 7331 much lower than in NGC 6946. Applying again the 
Kennicutt (1983) relation to the observed integral Hq flux (Young 
et al. 1996) after correcting this for intemal extinction in the way 
suggested by Devereux & Young (1993) we estimate a global star 
formation rate of 7 Mo/yr. This is about half the value found for 
NGC 6946. Furthermore, the gaseous and, in particular, the stellar 
disc (Broeils 1995) are more massive than in NGC 6946 so that the 
relative disc mass increase is only moderate. 

Finally we note that the adopted distances to the galaxies af- 
fect the dynamical stability condition of their gas discs (cf. Zasov 
& Bizyaev 1996). The surface densities and amplitudes of the ro- 
tation curves are distance independent. The epicyclic frequency 
scales inverse proportional to the distance, while the estimate of 
the turbulent velocity dispersion of the gas used throughout this 
study is distance independent. Thus assuming half of the presently 
adopted distance to NGC 6946, as suggested by de Vaucouleurs 
(1979), raises the stability parameter of the gas disc close to the 
threshold of dynamical instability. 



6 CONCLUSIONS 

We have formulated a local stability criterion of Galactic discs 
against axisymmetric density disturbances modelling the different 
dynamics of the stellar and gaseous components. The disc of the 
Milky Way is shown to be dynamically stable at all Galactic radii 
and probably over most of its past history. The inner parts of the 
disc of NGC 6949, a typical Sc galaxy from the Kennicutt (1989) 
sample, are found to be dynamically unstable. We have followed 
the ensuing dynamical evolution of the disc by numerical simu- 
lations. These show that such unstable discs evolve very rapidly. 
In order to stay in its present state the stellar disc would have to 
be effectively cooled by star formation. This seems to be actually 
observed in NGC 6946. However, the gas disc would have to be re- 
plenished by heavy accretion of gas, amounting to several times the 
present day gas disc mass during a Hubble time. 
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APPENDIX A: 

The stability of the self-gravitating composite gas and stellar disc is examined. Both disc components are approximated as infinitesimal thin 
sheets. In order to test the stability of the disc it is subjected to density perturbations of the form 

exp i{ujt + kR) , (Al) 

where uj and k denote the -complex- frequency and radial wave number, respectively. Since the most unstable wavelengths A = 27r/fc 
turn out to be small compared to the radial extent of the disc it is sufficient to study the stability of the disc in a localised theory (Toomre 
1964). Considering now a strip of the galactic disc, radial variations of the surface density and of the velocity distribution of the stars and the 
interstellar gas can be neglected. 



Al Stellar disc 

In order to describe the dynamics of the stars in the circular strip around the galactic centre we consider the collision-less Boltzmann equation. 



(A2) 



where / denotes the distribution function of the stars in phase space and H is the Hamiltonian. The distribution function of the un- 
perturbed disc, /o, is chosen according to equation (I) and normalised to the surface density Eq. Radial variations of fo are neglected. The 
- axisymmetric - disc response /i to a small axisymmetric potential perturbation, 5$ — <l>fc exp ifc/J, is calculated from the linearised 
Boltzmann equation 



21 

dt 



i + [/i,iJo] + [/o,5$] =0. 



In the following the plane stellar orbits are described according to the epicyclic approximation. The Hamiltonian Hq is then given by 



Ji 



A 



Ho = -R' + -Rtie - fio) - 2Ano{R - Ro), or alternatively Ho = ^ + —[J2 - QoRo) ~ Ti^oRo 



2 2 ' ^ ' ' ' 2 

where A and B denote Oort's constants. This leads to the equations of motion 



Ro {9 - riot) = m2 TTTT^ cos Wi , 



(A3) 



(A4) 



(A5) 



U = R = V2iZh 



COSUIi, 



V = Ro{e - Qo) + 2A{R - Ro) = 2B 



sm wi , 



(A6) 



where R, 9 denote polar coordinates. Ro is the mean galactocentric radius of the strip of the disc under consideration. is the mean 
angular velocity of the stars in the strip around the galactic centre. Ji and J2 are integrals of motion, Ji = (f/^ + ■f^V'^Y^'^ the radial 
action integral and Ro ( J2 ~Q,oRo) the angular momentum of a star referred to the mean radius of the strip, lui and 1112 denote the conjugate 
variables derived from the Hamiltonian Ho, 



A 



■W2 = — (J2 — QoRo)t. 
B 



Wl = Kt 

The Poisson brackets in equation (A3) are then given by 



[/o,(5$] = "TTT • = 5$lfccOSTOi. 

oJi OWi au 



(A7) 



(A8) 



Since the disc response is axisymmetric no df\/dw2 derivatives appear in equation (A^ ). This leads to the final form of the linearised 
Boltzmann equation 



dt dwi 



K 3\/3 
—2B 2TTafj 



v2«:Ji I iK<I>fc cos TOi exp ifc h 

<Ju J \ —2B 



sm Wl 



(A9) 



Assuming a time dependence according to equation (Al), 
/i, ^fe oc expit^;*, 



(AlO) 



equation ( A9 ) can solved by combining solutions of the homogeneous part of the equation and a particular solution of the inhomoge- 



neous equation found by 'variation of the constant', giving finally 



\/2kJ\ Oy 



2 sin (ttiaj/k) 



dw\ exp i 



\ ——k (sinuii + sin {wi + w'l)) 



(All) 
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Figure Al. 

curve. 



crit 

Curve of neutrally stable axisymmetric perturbations of a stellar disc with an exponential velocity distribution. Unstable solutions lie below the 



In order to have a distribution function /i, which is uniquely defined in velocity space, the integration constant of the solution (411) 
has been chosen in such a way that /i is periodic with respect to the angular variable wi (cf. Toomre 1964). The density response of the disc 
to the potential perturbation is found by integrating the distribution function /i over velocity space, 



Si. 



Jo Jo 2710-;^ \au J ^2kJx 



1 - 



2 sin {nuj/K.) 



dw'i exp i 



-Wi 



-k (sinuii + sin {wi + w'l)^ 



(A12) 



The surface density perturbations are assumed to be self consistent. Therefore the density and potential perturbations have to satisfy the 
Poisson equation, 

= 4ttGY;i5{z), (A13) 

where 5{z) denotes a delta functi on wi th respect to the vertical spatial coordinate z. If the wave numbers |fc] are sufficiently large, 2> 1, 
the Laplace operator in equation (A13) has nearly cartesian form and one obtains 



27rG|fc| 



(A14) 



Thus, the amplitude of the density perturbations, Ei^, cancels out of equation (A12). Finally one obtains after evaluating the quadratures 
with respect to Ji and wi the dispersion relation. 




UJ / K 



2 sm{-Kiu / k) 



dw'i 



oxp{i — w'i) 



l + |e(^)^cos2(^). 



(A15) 



where A^rit = 47r^GEo/«:^ is the critical wavelength and ^ = ^Tr'^aff / X^^.^^k'^ . For lo = equation (A15) gives the curve of neutrally 
stable perturbations in the (A, ^) parameter space as shown in Fig.|Al|. Exponentially unstable perturbations lie below the curve, whereas 
oscillatory solutions lie above. As can be seen from Fig. Al discs with ^ > 0.394 are stable at all wavelengths. This corresponds to the 
Toomre parameter given equation (2). 

In case of a three-dimensional disc of finite thickness the vertical oscillations of the stars have to be taken into account. In the epicyclic 
approximation, however, the vertical oscillations separate from the planar motions and the energy associated with the vertical oscillations has 
simply to be added to the Hamiltonian (A4). Vandervoort (1970a,b) has developed methods to solve the Boltzmann and Poisson equations 
for such a system. A rough estimate of the finite thickness effect on the stability of the disc can be found, however, according to Toomre 
(1964). It can be shown that potential perturbations of a disc with an effective scale height zo are reduced by a factor of 

1 — exp(— fczo) 



kzo 



(A16) 



with respect to the potential perturbations of an infinitesimal thin disc with the same surface density. In order to correct the dispersion relation 
for this effect, the right hand side of equation (A15) has to be multiplied by the reduction factor (A16). 
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A2 Gaseous disc 

The gaseous component of the galactic disc is modelled by an isothermal gas. As is well known (Binney & Tremaine 1987), the density 
response of the gaseous disc is given by 

= S^.o$fc^2_e2fc2+4f^„B' (A17) 

where Ego denotes the unperturbed surface density of the gaseous disc and c is the turbulent velocity dispersion. 
Inserting into equation (A17) the analogue to equation (A14) leads to the dispersion relation of axisymmetric perturbations of a gas disc, 

^ Acrit 1 ■ 1 ^ / A 1 n\ 

l = T-rT 7^^^ ' Q= (A18) 

|A| 1 (")^ I ^ ( ^cr.f -)2 ^ ttGEo 

The effect of the finite thickness of the disc can be taken into account analogous to equation (A15) by multiplication with the reduction factor 
(A16). 



A3 Stellar and gas disc combined 



In the case of a two-component disc the density perturbation in equation (A14) refers to the total surface density of the disc response. Thus, 
the density response of the gaseous disc (A17) has to be added to the density response of the stellar disc (A 12) and the sum then inserted into 
equation (A14). The resulting dispersion relation has the form 



1 = i^i- 1 f , '^^p^^-"'-^ I ':rf K (A19) 



where the parameters Acrit*, C* ™d ^crit.g, Qg are defined as in sections Al and A2 for the stellar or the gaseous discs, respectively. 
Setting Lu = gives the curve of neutral stable perturbations. 



© 0000 RAS, MNRAS 000, 000-000 



